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ABSTRACT 



The Hanle effect has been proposed as a new diagnostic of circumstellar magnetic fields for 
early-type stars, for which it is sensitive to field strengths in the 1-300 G range. In this paper 
we compute the polarized P-Cygni line profiles that result from the Hanle effect. For modeling 
the polarization, we employ a variant of the "last scattering approximation" . For cases in which 
the Sobolev optical depths are greater than unity, the emergent line intensity is assumed to be 
unpolarized; while for smaller optical depths, the Stokes source functions for the Hanle effect with 
optically thin line scattering are used. For a typical P Cygni line, the polarized emission forms 
in the outer wind, because the Sobolev optical depth is large at the inner wind. For low surface 
field strengths, weak P Cygni lines are needed to measure the circumstellar field. For high values 
of the surface fields, both the Zeeman and Hanle diagnostics can be used, with the Zeeman effect 
probing the photospheric magnetic fields, and the Hanle effect measuring the magnetic field in 
the wind flow. Polarized line profiles are calculated for a self-consistent structure of the fiow 
and the magnetic geometry based on the WCFields model, which is applicable to slowly rotating 
stellar winds with magnetic fields drawn out by the gas fiow. For surface fields of a few hundred 
Gauss, we find that the Hanle effect can produce line polarizations in the range of a few tenths 
of a percent up to about 2 percent. 

Subject headings: Polarization (Hanle Effect) - Stars: Magnetic Fields - Stars: Winds - Tech- 
niques: Polarimetric 



1. INTRODUCTION 

The Hanle effect is a diagnostic that refers to the modification of resonance line scattering polarization 
in the presence of a magnetic field. The effect begins to have an infiuence on the polarization at fairly small 
field strengths of just a few Gauss with sensitivity to fields up to around 300 Gauss. Experiments to describe 
the polarization from resonance line scattering date back primarily back to the first third of the 20th century 
(see Mitchell & Zemansky 1934). The influence of a magnetic field on the line polarization was explained 
first by a young physicist named Wilhelm Hanle. Hanle (1924) described the change of linear polarization 
by the magnetic field in semi-classical terms as arising from the precession of an atomic, damped, harmonic 
oscillator. From a quantum mechanical point-of-view, the effect is understood in terms of interferences that 
occur when the degeneracy of the magnetic sublevels in the excited state is partially lifted. The effect has 
come to have applications to many topics in atomic physics as described in Moruzzi & Strumia (1991). 

Astrophysically, the Hanle effect has been used only in the Sun. Applications include magnetic field 
measurements in the chromosphere, corona, and particularly in prominences and filaments (e.g., see Lin, 
Penn, & Kuhn 1998). A detailed description of the physics of the Hanle effect and polarized radiation 
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transport for solar studies appears in Stenflo (1994). However, for stellar astronomy, the Hanle effect remains 
relatively unknown. 

This paper is the fourth in a series to explore applications of the effect for other stars, especially hot stars 
with winds (Ignace, Nordsieck, & Cassinelli 1997 [Paper 1]; Ignace, Cassinelli, & Nordsieck 1999 [Paper II]; 
Ignace 2001 [Paper III]). It has been our goal to develop model line profiles with the Hanle effect that can 
be used to interpret spcctropolarimetric data for inferring the properties of circumstellar magnetic fields. 
The Zecman effect has had some success in measuring photospheric surface fields for a few hot stars that are 
not part of the extreme Bp class (e.g., Henrichs 2003; Donati et al. 2001, 2002; and Neiner 2002). However, 
Ignace & Gayley (2003) have shown that the Zeeman effect suffers serious technical challenges for use in 
diagnosing circumstellar magnetic fields from wind emission lines. For example, typical circular polarizations 
for such lines have peak values of only 0.01% for a field with a surface value of around 100 G. Fields of this 
magnitude can have a significant influence on the wind flow from hot stars (e.g., Maheswaran & Cassinelli 
1992; Babel & Montmerle 1997ab; CassineUi et al. 2002; ud-Doula & Owocki 2002). Small polarizations 
result for the Zeeman effect because the Zecman-split components that produce the circular polarization are 
incoherent and oppositely signed. The wavelength separation of the Zeeman split lines is small compared with 
typical astrophysical broadening effects. Consequently, the net circular polarization of the line is severely 
reduced owing to polarimetric cancellation from blending of the Zeeman components. The Hanle effect is 
a modification of the linear polarization that arises from resonance line scattering that does not suffer from 
polarimetric cancellation due to line blending. 

In a semi-classical description of the Hanle effect, the line scattering can be split into two parts: an 
isotropic scattering component and a dipole scattering component like free electrons (e.g., Hamilton 1947; 
Chandrasekhar 1960). The fraction of the scattering that is dipole- like is referred to as the "polarizability" , 
El, the isotropic scattering contribution is given fractionally as (1 — £^1). The introduction of a magnetic 
field induces a precession of the classical oscillator motion of the atomic scatterers. At the same time, the 
oscillations are damped, emitting their absorbed radiation at a rate given by the Einstein A-value for the 
transition of interest. Thus there is a competition between the radiative rate versus the Larmor precessional 
rate that has the frequency, ojl = eB/rrieC. Using these same arguments, Hanle (1924) deduced that a 
significant effect results when these two rates are somewhat comparable. We thus define a "Hanle ratio" 
given by 



B 2gL^L (B/b G) 

= .9i /-ins „-n ' (1) 
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where is the Lande factor, and Aui is the Einstein ^-value from upper level u to lower level I. The Hanle 
field strength Bnan is a conveniently defined parameter given by 



SHan = (2) 

2gLe 

that characterizes the magnetic field scale at which the the Hanle effect is important. Values of Bnan for 
typical strong resonance lines common to astrophysics are displayed in Figure 1. Note that the lines at 
shorter wavelengths tend to have larger Hanle fields, because A- values increase roughly as A~^. 

Recognizing that the Hanle effect is sensitive to modest field values of 10-100 G, we explored in Paper I 
its use for determining stellar wind magnetic field values and geometries from the total line polarization of 
optically thin emission lines. In Paper II we considered the Hanle effect in optically thin and resolved line 
profiles that from in equatorial disks. In Paper III expressions were determined for the depolarizing effect of 
having stellar radiation arising from a finite disk instead of from a point star. For thin lines the correction 
was found to be the same as that for electron scattering polarization as derived by Cassinelli, Nordsieck, & 
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Murison (1987). It was pointed out by Cassinelli, Nordsieck, & Ignace (2001) that in contrast with the solar 
case, the appHcation of the Hanle effect to hot stars is somewhat simpler and less ambiguous. This is for two 
reasons, (a) The hot star wind lines are highly NLTE so that virtually all bound electrons of metal species 
are in the ground state, whereas detailed NLTE considerations are required in solar applications, and (b) 
the stellar continuum radiation that is being scattered is relatively "flat" over the frequency width of any 
given line, while solar applications often involve scattering of chromospheric emission lines for which the flux 
changes strongly with wavelength across the resonance line profiles. 

The goal of these papers has been to move from an introduction of the Hanle effect for the broader 
stellar community to an investigation of its use by way of model calculations of increasing physical realism 

and numerical sophistication. In this paper we treat more realistic and more dynamically self-consistent 
magnetic field distributions in the wind, and we present results for the distribution of polarization across 
P Cygni line profiles. 

Section 2 outlines the line polarization calculations, and introduces a "single scattering" approximation 
that is used together with traditional Sobolev line profile analyses. We find that line optical depth effects 
have a significant influence on the polarized line profiles. Heuristic results are described in §3, including 
simple spherical shells with axial fields, and spherical winds with dipole fields. Results for a more realistic 
scenario for the field and flow are presented §4. The model chosen for this is the slow magnetic rotator 
field configuration based on the WCFields theory by Ignace, Cassinelli, & Bjorkman (1998). Significant 
line polarizations of around 1% are foimd to result. Lastly, §5 offers a discussion of our modeling, with 
implications for observations of hot star winds at ultraviolet wavelengths and future directions for modeling 
efforts. 



2. MODEL LINE PROFILES 

In the previous papers, we made calculations of the Hanle effect in emission lines that arc valid only in 
the optically thin case. Here we extend the results to allow for optical depth effects, by implementing escape 
probability theory and introducing the "single scattering approximation" for the emergent polarization from a 
volume sector. In this approximation the radiation at large optical depth is assumed completely unpolarized, 
and magnetically sensitive line polarization effects arise only in optically thin regions where a photon is 
considered to scatter only once. We make use of the Sobolev theory for resonance line scattering in stellar 
winds as developed by Sobolev (1960), Castor (1970), and Mihalas (1978). In Sobolev theory the radiative 
transfer is approximated as being local, because of the fast flow speeds and large velocity gradients that are 
present in stellar winds. The scattering and production of photons at any small volume in the medium can 
be assumed to have no interaction with any other portion of the flow. 

2.1. Standard Sobolev Theory 

In the standard Sobolev approach (e.g., as summarized in Mihalas 1978), the observable emission line 
profile is constructed from a ray-by-ray determination of the emergent line intensity for a fixed frequency 
within the line profile. An integration of these intensities is used to obtain the total emission or absorption 
at a particular frequency within the profile. The key variables for computing the line emission for a spherical 
flow are the velocity shift along the line-of-sight axis z, 



= -v{r) /U, 

where v is the spherical velocity flow, r is the radial distance from the star, and /x = 



(3) 

cosx is the polar 
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spherical angle from the observer's line-of-sight; the Sobolev optical depth 



where ki is the frequency-integrated line opacity, p = p{r) is the density, Aq is the line central wavelength, 
and a is related to the line-of-sight velocity gradient (as will be given in eq. [13]); and the emergent intensity 
along the ray at impact parameter p, 



IAp) = SAr) (1 - e-^^) , (5) 

where Si,{r) is the line source function. 

Because the flow speed is much greater than the thermal broadening only a small distance along the 
line-of-sight affects the intensity, and the integration over impact parameters at a fixed frequency in the line 
profile corresponds to an integral over an "isovclocity surface" . Such a surface is defined by v.^ = constant. 
A simplification implicit in equation (5) is that the surfaces of constant velocity-shift are single-valued with 
respect to intercepting rays of fixed p. More complicated cases involving non-local coupling associated 
with doublets or more general flow patterns are presented by Rybicki & Hummer (1978). Also common in 
applications of Sobolev theory is the distinct "core-plus-halo" approximation in which there is a contimiTim- 
producing photosphere and a distinct overlying outflow, in which the absorption of the continuum is computed 
via I*e~'^^ along each ray intercepting the photosphere. 

A different but equivalent approach based on escape probabilities was developed by Rybicki & Hummer 
(1983). The method involves integrating the emissivity for the line radiation over the emitting volume ac- 
counting for optical depth effects via escape probabilities. This approach is more suitable for our purposes 
not because of non-local coupling or complicated flow patterns (for we assume a radial velocity flow), but 
because we want to allow for axisymmetric magnetic field topologies viewed from arbitrary viewing incli- 
nations. It is convenient to solve for the emergent polarized flux based on a volume integration over the 
scattering medium. 

To describe our calculations, we begin with the basic escape probability approach to Sobolev theory. 
From Rybicki & Hummer (1983), the emergent luminosity from a volume element is 



dLi, = ^irji, pPfjt dr da dp,, (6) 

where jvp is the emissivity as given by 



jv p= KlpS^'ll){v - Ur,) (7) 

and Pij, is the directional escape probability given by 

= (8) 

for Ts = Ts{r,ii) the Sobolev optical depth as given previously. In equation (7), the factor if) is the line 
profile function. In Sobolev theory this function is traditionally a delta function; however, it is key only that 
this function be much more narrow than the wind-broadening of the line. For the fast winds of hot stars, 
this turns out to be an excellent approximation, since in the co- moving frame of the gas, V has a velocity 
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width of order v^h for a given ion species, and this is much smaller than the wind broadening that is of order 
Vao- The frequency appearing in the argument of ip is related to the velocity shift of equation (3) by 



1^0 



for the central wavelength of the line under consideration. The narrow profile response function ip 
effectively picks out, for a given sightline, a particular volume element or "cell" corresponding to that cell's 
Doppler shift with v = i/^. Hence the strong velocity gradients in hot star winds ensures that the radiative 
transfer for all volume elements is "local" , meaning all the cells are radiatively de-coupled from each other. 

The profile function ijj is defined such that 

Tpiu - u') du = 1. (10) 



f 

Jo 



For the purposes of numerical calculation, we consider a wind emission line profile to have a number of 
fixed frequency bins, each of width Au^ = (vo/cjAvz = Avz/Xq and in total spanning the full line width 
Avoo = (2t.'oo/c) i/Q- Then for a given volumetric cell existing in a spherical shell at fixed radius r, with fixed 
radial extent Ar, the luminosity contribution to the observed line emission will be 

AL^ = A-KKipS^Pu.r'^ArAa^. (11) 

Av^ 

Using the expression for the escape probability equation (8) and the Sobolev optical depth equation (4), and 
the relation between Au^ and Aw^, the luminosity from this cell becomes 

AL^ = 47r5^ (1-e-^s) (1 + a ijl^) vir) r Ar Aa ^ , (12) 

Auz 

where 

d\n.v ^ 

<j=- 1. 13 

ctlnr 

It is the form of the Sobolev approximation given in equation (12) that we use for computing emission line 
profiles. In particular, pure scattering resonance lines are assumed, for which the source function is 

Su = ^I., (14) 

Pesc 

where /? represents the escape probability averaged over the relevant solid angle. 

Pill) = ^ J^pi^dn = ^ J Pi^dfi. (15) 



In particular, the two (3's that occur in equation (14) are: the penetration factor 

r+i 



1 
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for fj,c = — Rl/r"^, and the escape factor 




+1 



(17) 



1 



It is useful to compare the thin and thick limits for the source function. For regions where rs <S 1 for 
all /U, one has from equation (8) that ^ 1 — 0{ts)- In this case /3c ~ 0.5(1 — iic) and /3esc ~ 1- Thus the 
ratio of /3's yields the well-known dilution factor 



The thick limit with ts ^ 1 is quite different. The escape probability becomes ~ Tg^ = f{r) [1 + 

(T(r)/i^], for /(r) a function of radius only. For a spherical shell of radius r, this escape probability factor 
is a parabola in /x, and thus a parabola also in Doppler shift in contributing to the line profile emission. 
For typical wind velocity laws that increase monotonically from some small value to an asymptotically large 
value, a will range from cr«— Itou^l. This implies that the parabolic factor can be either concave up 
or down, and this behavior affccs the line profile. The exact form of the source function thus depends on 
the velocity law v{r), and for our model line profile calculations, a typical beta wind velocity law, which to 
avoid confusion with the escape probabilities) will be written using 7, is adopted as given by 



where 7 is the velocity law exponent, and 6 < 1 sets the initial wind speed vq. For such a form, and with 
7 = 1, it can be shown that in the optically thick limit, the source function will vary as Si^ cx r~^. 

Having described the source functions for wind emission lines in the standard Sobolev approach, we 
turn now to a discussion of the polarizing properties of the Hanle effect. 



Computing the polarimetric properties of a beam of line radiation that is modified by the Hanle effect 
requires a number of coordinate transformations involving the local magnetic field at the scattering point, 
the star-centered frame, and the observer's coordinate frame. Figure 2 shows the stellar coordinates that are 
Cartesian (a;*, y*, 2:*) and spherical (r, ip). The solid point in Figure 2 represents a point of scattering, with 
the scattering geometry shown in Figure 3. There, the unit vector s for the scattered light is the direction 
toward the observer (i.e., along the z-axis). The spherical angle x is the angle of scattering from the radial 
direction, and a is an angle to the local field direction. Cylindrical coordinates {p, a, z) refer to the observer, 
with p being the impact parameter (not shown in the figure). We take the star to be axisymmetric, and 
tilted with respect to the observer's axis by an inclination angle i. 

The field topology is also assumed to be axisymmetric about the axis z«, and star-centered. In the star 
system, the orientation of the vector magnetic field at any point in the wind is given by the polar angles iJb 
and ipB- Polar angles 9s and Oi refer to the direction, measured from the local magnetic field axis B, of the 
scattered light and the radial direction, where the local radius vector f is taken to be the symmetry axis for 
the incident bundle of intensities from the stellar atmosphere. We will ignore limb darkening so that stellar 
intensities J* are the same from every point of the projected photosphere. Azimuthal angles 0s and are 
measured around B. It is only their difference that enters into the phase function that governs the Hanle 




(18) 




(19) 



2.2. The Hanle Source Functions for Single Scattering 
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effect, so we define 5 = (jig — '/'i, which is what appears in Figure 3. Details regarding the inter-relations of 
the various geometrical angles and their evaluation are given in Appendix A. 

The Stokes parameters /, Q, U, and V arc used to describe the polarization of the light, which will be 
represented as a "vector" , as for example the flux of light will be F = {Fj,Fq, Fu,Fv) = L/47r-D^. However, 
we shall always assume that V = (i.e., the circular polarization is zero), and so will ignore that component 
in what follows. 

For the Stokes source function, we have S = {Sj, Sq,Su). In Papers I-III, expressions were given for 
the source function elements, which we are repeated here in modified form: 



Si = Sii, 

Sq = Sqi cos 2% + Sm sin 2% , 
Su = -Sqi sm2is + Sui cos 2is , 



(20) 
(21) 
(22) 



where Sn, Sqi, and Sui are the source functions in the magnetic field system, and is is an azimuthal angle 
between that system and the observer's (see Fig. 2). The source function components are. 



Sn = J+-Ei {3K - J) 



Sqi 



:{1 — 3cos^ ^i) (1 — Scos^^s) +a{B,d) cosOg cos sin sin^j 



+ ^c{B,6) cos^^isin^^s 



- El {3K - J) 



^sin^ 0s (1 - Scos^ 6*0 



+a{B, 5) cos cos 0i sin 0s sin 9i — ■^'^{B, 5) sin^ 9i (1 — cos^ 9s) 



Sui = -Ei{3K-J) [-b{B, 6) cos 0i sin 0i sin 0s + d{B, 6) cos 9s sin^ 0i\ 



(23) 

(24) 
(25) 



where Ei is the polarizability, J and K are the standard Eddington moments of the radiation field (assumed 
to be determined by the incident stellar radiation field), and the local field strength enters through the 
functions a, b, c, and d, which are given by 



where the a's are defined by 



a = cosai cos(^ — ai), (26) 

b = cosai sin(5 — Q!i), (27) 

C = C0SQ!2 cos2(^ — 0:2), (28) 

d = cosQ!2sin2(5 — 02), (29) 



QLeB m B 

tanam = m = ttt; > (30) 

for the Lande factor, and m = 1 or 2. The case of non-magnetic scattering corresponds to = 0, for 
which the scattering source functions reduce to the expressions for pure Thomson (dipole) scattering except 
for El. In the strong or "saturated" limit of the Hanle efii'ect, « 7r/2, and soa = b = c = d = 0. 
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Finally, it is worth noting the analytic expressions for J and K in the current approximations. Assuming 
a star of uniform brightness, the Eddington moments are given by 

J = i- y I^dn = W{r)h, (31) 

and 

K = ^ j hn''dn = W{r) h (^l±i|±ifl^ . (32) 

The factor ^(3if — J) is the familiar finite star depolarization factor of Cassinelli et al. (1987). Prom the 
preceding two expressions, this geometric factor is given by 

i(3i^- J) = W^(r)4 X (l + /xe). (33) 

At large radius the factor scales like r~^, which is the point source approximation. However, right at the 

photosphere, the it drops to zero. This results because the symmetry of the radiation field over half the sky 
as seen by a scattering particle yields no net polarization of the scattered light. On the other hand, limb 
darkening prevents the factor from entirely vanishing at the photospheric level. Although we will ignore 
limb-darkening in our models, the effect is relevant for the use of the Hanle effect in studies of the sun that 
probe scattering geometries around the solar photospheric level (e.g., see Stenflo 1994). 



2.3. The "Single Scattering" Approximation for Polarized Emission Lines 

It is well-known that multiple scattering of linearly polarized light generally leads to a depolarization. 
The polarized intensity from most media of astrophysical interest is small, and so one can approximate the ra- 
diative transfer as predominantly unpolarized. The depolarization that arises from multiple scattering would 
suggest that an emergent intensity beam derives the bulk of its polarization from its last scattering event. 
In the "last scattering" approximation, radiation propagates through a medium as a beam of unpolarized 
radiation until it emerges, with a sense of polarization appropriate for its last scattering interaction. 

For the wind case, it is a "single scattering" approximation that is the picture most relevant for our 
analysis. In Sobolev theory the radiation transport for the line emission is localized to relatively small 
volumes, sometimes referred to as "resonance zones" . In the terminology of Monte Carlo radiative transfer, 
one can envision many photon "packets" of stellar light that penetrate into these zones, scatter several 
times, and emerge toward an observer. Owing to the probabilistic nature of treating the scattering, the 
scattering history of each packet will be different. Although each will be characterized by a polarization, 
the polarization position angles will on average be random, thus summing to small net polarizations. In 
contrast photon packets will typically scatter only once inside optically thin resonance zones. These will 
all have a polarization, but more importantly they will on average all have the same polarization position 
angle so as to sum constructively. Consequently, we employ the simplifying assumption that the scattered 
light from optically thick (rg > 1) resonance zones is completely unpolarized, whereas light emerging from 
optically thin regions (rg < 1) is treated as if single-scattered, for which we use the Hanle source functions 
as previously defined. 
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3. RESULTS FOR AXIAL AND DIPOLE MAGNETIC FIELDS 

In this section we present polarized profile results for heuristic purposes, focussing on two simple mag- 
netic geometries: an axial magnetic field and a dipole magnetic field. In each case the underlying wind is 
assumed to be spherical. This latter assumption is not dynamically consistent with the chosen field topolo- 
gies, but we adopt the condition so as to isolate the influence of the Hanle effect for the polarization of wind 
emission line profiles. 

The viewing perspective is described completely by the inclination angle, i in the stellar frame. Since 
many of the strong resonance lines common to astrophysics are Li-like doublets, we use Ei = 0.5 for our 
calculations, which is the value appropriate for the shorter wavelength component for many of these doublets 
(e.g., Chandrasekhar 1960). 

The line strength will be characterized with an optical depth scale given by 



PO ■^O /q.N 



and so the Sobolev optical depth can be expressed as 



Ts=Tiw-^x-^{l + fi^a)-\ (35) 

with w = v/voo the normalized radial velocity and x = r/R>, the normalized radius. Some example P Cygni 
profiles as calculated with our numerical code are shown in Figure 4 for a range of optical depth scales r; 

and with a 7 = 1 velocity law. All of our model profile calculations will involve 7 = 1, and the case of 
Ti = 1 will typically be adopted. For simplicity we further assume the ionization of the wind is constant with 
radius. The ionization of different atomic species is generally a function of radius (e.g.. Drew 1989); however, 
some ions can be; dominant over a substantial range of radii. Moreover, although variations in ionization can 
change the level of polarized fiux, the percent polarization of the profile is not altered when ionization varies 
with radius only. For this reason we concentrate on the constant ionization case to emphasis the influence 
of the Hanle effect for the polarized line emission. 

The different magnetic wind models will be parametrized in terms of the ratio of the surface field 
strength relative to the Hanle field -Bnan- Since the figures will display polarized line profiles for a range 
of ratios B^/Buan, one may view the results in two different fashions. For a line with a given Hanle field, 
different polarized profiles represent how the line would appear for stars of different surface field strengths. 
Conversely, for a star of a given field strength, the different polarized profiles represent lines of different 
Hanle fields from the same wind (i.e., assuming that all the lines form over the same spatial region). It is 
useful to scrutinize the results that follow with both perspectives in mind. 



3.1. The Axial Magnetic Field Case 

For an axial magnetic field, there is a tremendous simplification of the Hanle effect problem in terms 
of geometry, because the magnetic field B = is everywhere parallel to the assumed axis of symmetry 

for the star with i)b = and the field strength is taken to be constant. This then is the easiest case to 
consider for the Hanle effect. In fact, the expressions describing the Stokes Q and U source functions for 
the light scattered by an optically thin and uniform ring of matter whose symmetry axis passes through the 
origin is analytic. For a spherically symmetric wind, such a ring represents the fundamental building block 
for constructing theoretical emission line profiles. 
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For the case at hand - spherical symmetry in the expansion and density - it is the scattered hght as 
integrated around the ring that is important for our consideration of unresolved sources. Consider a ring 
whose axis is directed toward the observer. All points on the ring have the same Doppler shift. Note that the 
axis of the ring is not the same as the symmetry axis describing the magnetic field. Using equations (23)- 
(25) with 9s = i, 6i = i), S = —(p and ig = 0, the azimuthally integrated (i.e., in a) Stokes Q and U source 
functions for axisymmetric rings are found to be 



Sq = --El sin^iSoir) (l-3cos^x) [(1 - 3cos^ i) + 4cos^ ai cos^ i - cos^ a2 (1 + cos^ i)] , (36) 
8 

and 

Su = ^Ei cos i sin^ i So{r) (l — 3 cos^ x) sin 2ai + ^ sin 2a2^ , (37) 

where So{r) = 3K — J is a function of radius only, and the Hanle mixing angles ai and a2 can vary with 
location (r, i?) . Assimiing for illustrative purposes that the field strength varies with radius only, several 
interesting points can be made about the Hanle effect in omission lines. 

Consider an optically thin spherical scattering shell. Each isovclocity contour is a ring centered on the 
observer's viewing axis and distinguished from other rings by the angle %. Figure 5 shows how the Q and U 
source functions will vary from ring to ring, with a constant viewing inclination for the curves at left, and a 
constant field strength for the curves at right. 

Even though both the ring and field are axisymmetric, note that a J7-signal persists. This is because of 
the field's handedness in producing the Hanle procession. Notable points include: (a) Both Q and U vanish 
for the case B = 0, since the non-magnetic scattering polarization vanishes for a spherically symmetric 
density distribution, (b) The polarized Q and U source functions both vanish at the "Van Vleck" angle, 
cos^X =1/3. (c) The i7-signal scales with inclination such that the [/-polarization vanishes for pole-on and 
edge-on viewings. (d) For spherical expansion the [/-profile integrates to zero across the line, (e) Finally, in 
the saturated limit of the Hanle effect, Su =0 and only a Q-signal remains. Indeed, in this limit the source 
function Sq is especially simple with the form 

Sq = ~Ei sin^ i 5o(r) (l - 3 cos^ x) (l - 3 cos^ i) . (38) 
8 

In addition to = for the two isovelocity rings situated at cos^ x = 1/3, the Q source function also 
vanishes for viewing inclinations at cos^i = 1/3. 

Figure 6 shows the polarized line profile (eq. [12]) for this geometry from our code, assuming a constant 
field strength and using the single scattering approximation with 7 = 1, r; = 1.0, and Ei = 1. Displayed are 

the continuum normalized Q and U Stokes flux profiles for an axial field and a spherical shell of fixed radius 
in percent polarization for the same range of viewing inclinations (relative to the axial field) as in Figure 5, 
but with a different selection of field strengths to better distinguish between the model line profiles. The 
vertical plot scale is different because the curves for Figure 5 are source functions, whereas those in Figure 6 
arc line fluxes. Since we explicitly assume that polarized line emission emerges only from optically thin 
portions of the wind, one might expect that the profiles that are evaluated numerically would resemble those 
of Figure 5, at least qualitatively. And indeed the numerical polarized line profile and the analytic source 
function do show some similarity, at line center, but whereas the analytic derivation predicts relatively strong 
polarization at the line wings, the polarized flux drops to zero for the numerical calculations. In particular, 
note that the line- integrated U fiux does not necessarily vanish, contrary to expectations. 



- 11 - 



These differences between the source function and line profile can be understood with the help of 
Figure 7, which shows an illustration with the star at center and dotted lines for isovelocity zones. The 
observer is along the z-axis located to the right. The solid curves are the locus of points with Sobolev optical 
depth rs = 1 (assuming a velocity law with 7 = 1 and line optical depth scales ti as indicated). Interior to 
sets of paired solid curves, the line optical depth is greater than unity, and exterior the line optical depth 
is smaller than unity. The heavy dark curves are for t; = 1, which represents the case most typical of 
our line profile calculations. In this case the figure indicates that there is no spherical shell that can be 
drawn which does not intercept the region of rg > 1. Moreover, the interception, for shells of large radius, 
occurs preferentially at those parts of the shell that arc primarily fore and aft of the star with respect to the 
observer, thus corresponding to Dopplcr shifts in the line wings. 

Importantly, we discover that there are two new effects that arise from a consideration of line optical 
depth effects. First, although spherical symmetry is used, the Sobolev optical depth for fixed radius varies 
in going from the front side of the shell to the rear side, with some portions being optically thick and others 
optically thin. Second, for some shells, every point can be optically thick. For example in the case r; = 1, 
Figure 8 shows that every point inside r « 2.5i?* has Tg > 1, and so the Hanle effect is effectively "blind" 
to the circumstellar magnetic fields between the wind base and this radius. 

This latter point implies two interesting corollaries. First, lines of different strength (as characterized 
by the different ti values in the figure) will allow the Hanle effect to probe magnetic fields to different radial 

proximities to the wind base. Second, the Hanle and Zceman effects can be used in tandem. For a star such 
as 9^ Ori C with a strong multi-kiloGauss surface magnetic field, the Zeeman effect as applied to photospheric 
lines yields information about the surface field distribution and strength. The Hanle effect as applied to wind 
lines will sample the circumstellar magnetic field, which will be weaker than the surface field. Together, the 
line polarizations might be used to reconstruct how the magnetic field and wind flow affect one another. 

In summary, considerations of the Hanle effect in a wind flow with an axial magnetic field reveals that 
substantial polarizations of up to 4% can result near line center. The Stokes Q and U profile shapes are 
characteristically symmetric, and the polarization position angle does not vary across the profile. A point 
that is significant for observational studies is that the line-integrated polarized flux docs not in general 
vanish, implying that a net polarization could be measured even for poorly resolved line profiles. The axial 
field topology is especially simple, and the derived peak polarizations are to be considered as the best possible 
because a constant vector magnetic field was assumed everywhere in the wind, whereas in real winds the 
field strength will generally diminish with radius and the field orientation will normally vary from point to 
point. The major result of this section is that using the axial field scenario as a control case, we find that 
line optical depth effects do substantially modify the shapes of the polarized profiles and the portions of the 
wind that can be probed by the Hanle effect. 

3.2. The Dipole Magnetic Field Case 

Now we apply our theory to extended and radially outflowing envelopes that are threaded with a dipole 
magnetic field. The combination of a dipole field and a radial wind are not compatible, since either a strong 
dipole field will result in wind confinement (Babel & Montmerle 1997ab; Cassinelli et al. 2002) or a strong 
wind fiow will drag out the magnetic field, distorting the dipole into a more nearly radial geometry (ud-Doula 
& Owocki 2002). Still the dipole case represents a more complicated field topology than offered by a simple 
axial field and so is worth investigating. Also, this is the topology used in section 3.2 of Paper II, that had 
a more simplistic assumption for the effects of optical depth. Here we see the influence of optical depth in a 
better treatment of the radiative transfer. 

It should be mentioned that a dipole magnetic field is purely radial at the poles, thus by symmetry 
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yielding no Hanle effect. In the equatorial plane, the field is purely axial; everywhere else, the field is a 
vector sum of radial and latitudinal components. At no latitude is there a toroidal component of the field. 
In the next section, the influence of a toroidal fleld will be explored in the context of winds from slowly 
rotating stars. 

Polarized line profiles have been calculated for a range of line strengths, viewing inclinations, and Hanle 

ratios, with results displayed in Figure 8. The two panels at left are for edge-on views with i = 90°. In this 
case Stokes U is identically zero throughout the line profile owing to the symmetry, and so only Stokes Q 
is shown (normalized to the continuum flux and plotted in percent). The profile polarization P and the 
variation of polarization position angle 9p are shown in the right panels. The Hanle ratios quoted in the 
figure refer to i?,/i?Han; however, we remind that the field strength of a dipole diminishes with radius as 
B oc r^"^, and so a range of radii and field strengths contribute to the observed line polarization. 

Rciturning to panels at left, the upper one shows that stronger Hanle ratios lead to lines of higher 
polarization. The left-right symmetry in the polarized profile corresponds to a fore-aft spatial symmetry 
with respect to the star. The panel below illustrates the effect of line optical depth, with solid for r/ = 0.1, 
short dash for r; = 0.3, long dash for ti — 1, dash dot for t; = 3, and dotted for r; = 10. For small optical 
depths, the polarized line emission arises from deep in the wind. As t; is made to increase, the region of 
> 1 moves outward, and the polarization peaks drift toward larger Doppler shifts in the line. 

Initially, increasing r; leads to larger line polarizations, because there are more scatterers. For larger 
optical depths, the region of > 1 has swelled sufficiently in extent that the overall flux of polarized line 
emission begins to drop. Interestingly, at low values of t;, stellar occultation results in an asymmetric profile 
of polarized line flux, with the redshifted emission somewhat suppressed. For larger t;, the proflle becomes 
symmetric. This can be understood in terms of Figure 7. For large t; , all the points lying directly behind 
and in front of the star do not contribute any polarized line emission, because those points all have rg > 1, 
hence occultation has no influence; at small r; this is not case, and so occultation makes a difference. 

Returning to Figure 8, the two panels at right show the effects of viewing inclination, which in the case 
of a dipole field, yields both Q and U profiles. The inclination values are indicated. The upper plot is the 



total polarization with P = yPg + Fu/ -^tot , and the lower plot is for the polarization position angle defined 
by tan2^p = Fu/Fq. Note that Ftot includes the redshifted P Cygni line emission and the blueshifted 
absorption, making the polarized profile in P asymmetric. On the whole the line polarization is seen to drop 
with viewing inclination, similar to the sin^ i rule of Brown & McLean (1977) for axisymmetric Thomson 
scattering envelopes. 

The variation in the polarization position angle is especially interesting. Our assumption is that the 
polarization arises from the optically thin portions of the wind. This suggests that Fq and Fu are roughly 
proportional in terms of optical depth effects. So the ratio Fu/Fq that determines the polarization position 
angle 6p is to first order independent of optical depth. Its amplitude is set by the viewing inclination, and 
its variation with Doppler shift is governed by the Hanle effect. The line types for 9p in the lower right plot 
of Figure 8 are the same as in the panel above it. In the absence of a magnetic field, the spherical symmetry 
ensures that both Q and U are zero, so that it is the Hanle effect that produces the line polarization and the 
position angle variations both. If the envelope and velocity field were instead axisymmetric, then position 
angle rotations could result even without a magnetic field. However, the Hanle effect would distinguish itself 
by virtue of how the position angle variations differed between lines of different A- values and optical depths. 

In comparison to the polarized profiles from spherical winds that were modelled in Paper II (Fig. 8 

in that work), the profile results presented here are seen to be qualitatively similar but to differ somewhat 
quantitatively. Under the single scattering approximation, the polarized profiles have significant polarization 
near line center, whereas those of Paper II were more centrally depressed. Although the polarizations in 
Paper II were overestimated by about a factor of 2, the overall asymmetric profile shape in P remains because 




- 13 - 



of the normalization by the P Cygni profile. 

To sum up, the main lessons to draw from the dipole magnetic field case are that overall stronger line 
polarizations result from larger Hanle ratios, the strongest line polarizations are for lines that are relatively 
thin, the polarization is reduced for viewing inclinations that are further from edge-on, and the influence 
of the Hanle effect is betrayed through the observation of position angle rotation effects in the line profile. 
Having established some insight into the Hanle effect of magnetized winds, and a general sense of the line 
polarizations that can result (up to about 0.7%), we next consider a more realistic case in the context of 
winds from slow rotators that drag out a relatively weak surface magnetic field. 

4. RESULTS FOR A SLOW MAGNETIC ROTATOR 

In this section, in order to define a realistic but tractable magnetic geometry, we assume that in the 
case of weak magnetic fields, the star will have a "slow magnetic rotator wind" . This means that the field 
does not play a role in accelerating the flow to terminal velocity. We are motivated by the fact that the 
Sun is a slow magnetic rotator (Belcher &: MacGregor 1976). The overall weak solar magnetic field does not 
dominate the wind flow, but is believed to have contributed to a long-term braking of the Sun's rotation (as 
first analyzed by Weber & Davis 1967). The field produces some transfer of angular momentum to the wind, 
and in the process leads to some rotational distortion of the flow. For simplicity we also assume that the 
rotation is slow enough that the wind velocity field is basically spherical, and that the field is dragged out 
into "streak lines" by the flow. Although the velocity field is treated as a function of radius only, we shall 
consider both spherical and axisymmetric density distributions for the line polarization and the influence of 
the magnetic fields. 

For the underlying field topology, we make use of the WCFields model (Ignace et al. 1998) that describes 
an initially radial field that is dragged out by the wind from a rotating star. The model is based on the Wind 
Compression theory of Bjorkman & Cassinelli (1993). Expressions governing the wind flow and magnetic 
field geometry are summarized in Appendix B. The salient features are that (a) a given wind model is 
predominantly determined by the ratio of the equatorial rotation speed v^ot relative to the wind terminal 
speed Voo, (b) the field is initially assumed to be radial at the wind base (i.e., a split monopole), (c) the 
fields are "weak" and frozen-in so that the wind flow geometry determines the magnetic geometry, and (d) 
asymptotically, the field becomes dominantly toroidal at equatorial latitudes while remaining largely radial 
near the pole. 

As described in Appendix B, we adopt approximate expressions to describe the flow and magnetic field 
that can be used for slowly rotating stars (i.e., Wrot/^^oo !)• In this limit spherical geometry for both the 

wind velocity and density can be adopted . As an example, a slow rotator model that can produce interesting 
line polarizations is the case v^ot/voo = 0.08. For a terminal speed of 2000 km s~^, this would correspond 
to a rotation speed of 160 km s~^, which is not atypical of observed vsmi values in O stars (e.g., Penny 
1996), and amounts to rotating at about 25% of break-up. The asymptotic density contrast between the 
equator and pole, Poq/Ppoio, in this case is merely 1.25 (i.e., with 7 = 1). Initially, we shall ignore this density 
variation in our models but will return to the issue of how a wind density that deviates from spherical affects 
the line polarization. 

4.1. WCFields with a Spherical Wind Density 

Polarized flux profiles have been computed for the cases Vi-ot/voo = 0.03, 0.08, and 0.13, and the resiflts 
are plotted as percent polarization in Figure 9. Only the Q profiles are plotted. There is a [/ profile, but the 
level of polarization is about an order of magnitude less than in Q, and so we do not show it. In each panel 
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the different curves arc for different Hanle ratios, witli -B*/-BHan = 0, 3, 10, 30, and 100, such that stronger 

polarizations result for larger Hanlc ratios. 

For all of these cases, the polarized profiles are seen to be single-peaked, in distinction to the dipolc field 
case. The profiles also show a net negative polarization, implying a position angle that is orthogonal to the 
axis of rotation. This results from the significant toroidal magnetic field component, and can be understood 
as follows. At the base of the wind, the field is initially a split monopole, with no toroidal component. Both 
toroidal and latitudinal components begin to develop as the flow draws the field out in radius (see Fig. 10 and 
App. B). The radial component of the field decreases asymptotically as r~^, whereas the toroidal component 
will decrease as sin^?, and the latitudinal component as sin2i?, for i? the stellar co-latitude. At radii 
where the Sobolev optical depth is thin, the latitudinal field component is negligible (although it docs account 
for the small level of U polarization), and the toroidal component is starting to become comparable to Br. 
The field is mainly radial near the poles, so that in these regions the Hanle effect is minor. Around the 
equator the field is becoming more and more toroidal. 

Consider the isovelocity zone corresponding to line center, which is the plane of the sky as intercepting 
the center of the star. Without a magnetic field, the polarized intensity is symmetric about the linc-of-sight 
to the star. By our convention the polarization is negative around the polar limb, and positive around the 
equatorial limb. The field is radial near the pole, so the polarization of scattered line radiation does not 
change in that region. However near the equator, the toroidal field is into the plane on one side, out of it on 
the other. This leads to the geometry of the classic laboratory Hanle effect, and in the saturated limit, the 
scattered light becomes completely depolarized. Integrating around a circle that is centered on the star and 
lies in this plane, it is clear that the net polarization will no longer be zero, but negative, being dominated 
by the un-modified polarized fiux from the polar zones. Similar kinds of arguments can be used to show that 
the scattered light from all isovelocity zones will tend to be negative. 

Figure 9 shows that the line polarization is stronger for more rapidly rotating stars. This just reflects the 
fact that in WCFields theory, ~ I'rot/^^oo, so that for a given location in the wind, faster rotation implies 
that the ratio B^/B^ at a given radius is larger. Figure 10 shows the variation of the field components with 
radius for a stellar co-latitude of = 45°, chosen such that the peak value of B^ will be roughly maximized. 
The field components are plotted as normalized to the total field strength Btot = + + -B J . Since B^ 

dominates over most of the radii shown, the overall field strength drops roughly as r^^, although somewhat 
less rapidly at the larger radii owing to the increasing importance of B^p. As noted previously, at very large 
radius where the toroidal component dominates, the field will decrease only inversely with radius, instead 
of quadratically. This figure is for the case of J'roi/uoo = 0-08; for faster rotations both the latitudinal and 
toroidal components would be relatively larger, and vice versa. 

Especially important to realize is that the Hanle effect is a non- linear effect. Although the bulk of the 
scattered light that produces the polarized line fiux comes from radii where rg is of order unity and somewhat 
less. Figure 10 indicates that the field is predominantly radial at these locations which should yield no Hanle 
effect. Thus it is the non-radial components that lead to a Hanle effect in terms of introducing an asymmetry, 
and yet the modification of the scattering polarization is determined by the total field strength, and so the 
radial field is still relevant. Consequently, the Hanle effect is in principle sensitive to the full 3D magnetic 
field, in terms of field strength and field geometry. The main diagnostic challenge is the fact that the net 
line polarization at any frequency in the profile is a strongly convoluted property. At the very least, it is 
evident that larger Hanle ratios lead to greater polarization. Next we discuss observational prospects and 
diagnostic approaches. 

Using the case of t;rot = 0.08t;oo, we have examined the influence of line optical depth and viewing 
inclination in greater detail. Figure 11 shows a series of model line profiles at this fixed rotation speed 
as T; and i are varied. The uppermost panel shows the effect of line optical depth for the case i = 90°. 
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Interestingly, line optical depth had a stronger influence on the line polarization in the Dipole fleld case. 
Recall that as r; is made to increase, the region over which the line is optically thick grows in extent, and 
this region does not contribute to the line polarization. The field strength of a dipole decreases as r~^, such 
that for r; > 1, the Hanle eflfect samples the magnetic field where even at modest radii it is quite small, so 
that the polarization is also small. The WCFiclds case is different. The field strength drops initially as r^^ 
for a split monopole, slowly transitioning to a predominantly toroidal field that decreases only as r~^. This 
is much more gradual than for a Dipole field, so that the line polarization does not drop significantly until ti 
becomes exceedingly thick r; > 3. Indeed, the line polarization is larger for a somewhat thick line at r; = 10 
than for a thin line at rj = 0.3, because at low optical depth, most of the optically thin line radiation comes 
from the inner the wind, where the field is largely radial with almost no Hanle effect. 

Figure 11 also displays the effect of viewing inclination in the lower four panels: the two at left showing 
polarized flux in Q and U, and the two at right showing the total degree of polarization P and the polarization 
position angle 6p. In each panel the line type corresponds to the same inclination angle value as indicated 
in the panel for Fq. 

4.2. WCFields with an Axisymmetric Wind Density 

A net resonance line scattering polarization occurs even when there is no Hanle effect (i.e., B = 0), if 
the distribution of scatterers is asymmetric. The scattering is similar to Thomson scattering, as previously 
noted, reduced by the factor Ei, but increased by the much larger cross-section compared to Thomson 

scattering. 

The previous section ignored this polarization contribution by assuming the wind density to be spherical. 
The motivation was to isolate the influence of the Hanle effect. Now relaxing this approximation. Figure 12 
shows polarized line profile shapes for the case Urot/^^oo = 0.08 and 7 = 1, which has an equator to pole 
density contrast of 1.25. The upper dark line is the polarization from pure resonance scattering with no 
Hanle effect. The dotted curves are the profiles with the Hanle effect, to be compared with Figure 9 that 
assumed an underlying spherical wind density. So the density distribution does introduce a kind of bias 
for the line polarization for its interpretation in terms of the Hanle effect. However, it is not a large effect, 
and a strategy that targets lines with different Hanle fields should be able to correct for the influence of 
non-spherical density. 

Allowing for p(r, fi) also implies that scattering polarization by free electrons in the wind may complicate 
the interpretation of the line polarizations. Figure 13 shows the expected continuum polarization from mildly 
distorted stellar winds as viewed edge-on. For this figure we are using a sequence of Wind- Compressed Zone 
(WCZ) models (Ignace, Cassinelli, & Bjorkman 1996) in conjimction with the expressions of Brown & 
McLean (1977) for the polarization from optically thin axisymmetric envelopes. The results also include the 
finite star depolarization factor (Cassinelli et al. 1987). The lower axis is the ratio Urot/^^oo) and the upper 
axis is the asymptotic value of the cquator-to-pole density contrast. Electron scattering depth of tq = 0.03, 
0.1, and 0.3 for an eqiuvalent spherical wind were used as reference models. The three points indicate the 
rotation values used for the Hanle effect models presented in the previous section. Relatively few O stars 
are known to be intrinsically polarized at the 0.1% level, but typically tq < 0.1, so this is not inconsistent 
with the small distortions implied in our models (McDavid 2000). 

In any case continuum polarization of any kind, including interstellar polarization, is not a significant 
complication to our analysis: It turns out that for the strong Li-like doublets commonly observed in the 
UV and FUV band, the red component has Ei = 0, meaning that the line scatters isotropically, making 
no contribution to the line polarization and thus not susceptible to the Hanle effect. Any polarized flux 
observed in the red line of the doublet must therefore arise from continuum processes, and since the doublet 
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components arc close together in frequency, the polarization of the blue component which can show a Hanle 
effect may be straightforwardly corrected for the continuum polarization. 

4.3. Field Topology 

The Hanle effect is sensitive both to the field strength and the magnetic geometry. This sensitivity 
can be made especially graphic. Wc previously noted that for the WCFields models, there is a small U 
polarization. In Figure 14, wc plot the polarized profiles for a model with Wrot/i^oc = 0.08, i — 90°, and a 
spherical wind, but now for a magnetic field that is initially outward radial at the wind base (i.e., a magnetic 
monopole - unphysical, but used here only to make a point). The result is that the Q profiles are much 
reduced in scale, and a fairly substantial anti-symmetric U profile results. The models arc identical in all 
respects to the corresponding models of Figure 9, except that the radial field in the lower hemisphere is 
outwardly directed instead of inwardly directed. Thus we find that it is still the case that larger Hanle 
ratios yield stronger line polarizations, but moreover, that the U profile is sensitive to the overall symmetry 
(top-bottom and left-right) of the global magnetic geometry. 

In practice, the symmetry axis of the magnetic field (if one exists) will not be known a priori. The 
Stokes Q and U line profiles for a source arc measured according to observer-defined axes. It may be possible 
to identify symmetry in the magnetic field by rotating the observed polarizations using Mueller matrices to 
look for symmetry in the polarized profiles. Generally, a constraint on the symmetry of the circumstellar 
magnetic field would be estimated as part of the profile modeling. 

5. DISCUSSION 

This contribution has focussed on the Hanle effect in resonance scattering lines common to hot star 
winds. We have stressed the fact that the Hanle effect is non-linear, in that it is the non-radial field 
components that give rise a change in the line polarization, yet at the same time, it is the total ficild strength 
(including the radial component) that governs the amount of the change. Consequently, the Hanle effect can 
be used to infer the full 3D magnetic geometry in the wind. For a given set of scattering lines, stronger fields 
will lead to higher line polarizations. Field topology also has an influence, as evidenced by the double-peaked 
polarized profiles for a dipole magnetic field, versus the single-peaked profiles for slow magnetic rotators. 

The Far Ultraviolet Spectro-Polarimeter (FUSP) is a sounding rocket payload with the goal of mak- 
ing the first spcctro-polarimctric measurements in the wavelength range of 1050-1500 A for several stellar 
targets (Nordsieck 2003). With a 50 cm primary, the instrument will have a spectral resolution of 0.65 A, 
corresponding to a resolving power of A/AA = 1800 at a wavelength of 1170 A. This in turn corresponds to a 
velocity resolution of about 170 km s""'^, which is sufficient to resolve wind-broadened lines with typical full 
widths of a few thousand km s~^. Targeting the star C, Ori, which was observed in the FUV with Copernicus 
(Snow & Morton 1976), it is anticipated that FUSP will produce the first detection of the Hanle effect in a 
star other than the Sun. 

To correctly interpret the polarization of scattering lines that will be measurable with FUSP, care must 

be taken to identify all sources of polarized radiation. For example, interstellar polarization is endemic to 
all polarimetric studies of distant objects, and there are standard techniques for its correction, based on 
wavelength dependence and lack of variability. 

Extracting the influence of a non-spherical wind geometry for the line scattering polarization can be 
approached in two different ways, (a) If the interstellar polarization can be determined, the polarization 
outside line frequencies due to Thomson scattering can be used to model the underlying wind geometry. 
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However, this is difficult in part because it may not be clear that the line forms in the same region that 
gives rise to the electron scattering polarization, (b) Alternatively, since B(r) is fixed for any given star, one 
can use a multi-line approach to infer the polarization arising from an aspherical distribution of scatterers. 
Different lines have different Hanle field values, so that lines with large Hanle values can be used to estimate 
the non-magnetic contribution. Once again, care must be taken since different lines may not form in the 
same spatial locations; however, the P Cygni profiles themselves could be modelled to obtain information 
about the ionization distribution in radius and latitude. 

The study of the Hanle effect for circumstellar environs and the development of radiative transport 
techniques is ongoing, especially in terms of relaxing the single scattering approximation for the line polar- 
ization. In this regard one can use expressions in Jeffery (1990) for "Sobolcv-P" theory to show that the 
single scattering approximation captures the flavor of the polarimetric behavior of scattering lines, but that 
in fact there is non-zero polarization from line scattering even at rather high Sobolev optical depths of a 
few, and even at optical depths of just a few tenths, the line polarization is not truly described by single 
scattering. We anticipate using Monte Carlo radiation transport techniques in the future to model the line 
polarization more accurately. The proven versatility of these techniques will be useful in handling more 
realistic wind models. 

Another effect ignored here needs to be explored. While our analysis applies directly to resonance 
transitions in ions with no nuclear spin (e.g., O VI, Si IV, C IV, Mg II, and Ca II), in other Lithium- like ions 

such as P V, N V; Na I and K I, the presence of nuclear spin breaks the degeneracy of the spin states, and 
this alters the value of the polarizability Ei making it sensitive to optical pumping. In fact in the presence 
of optical pumping, the magnetic field has another effect on the line polarization through "magnetic re- 
alignment" (Nordsieck 2001), or in Solar Physics parlance, the "second Hanle effect". We are investigating 
the regime in which the Hanle effect and the magnetic re-alignment physics both apply. 

We thank Jon Bjorkman for useful discussions regarding line polarization effects. We are especially 
grateful to the referee, Marianne Faurobert-Scholl, for several helpful comments. Support for this research 
comes from a grant from the National Science Foundation (AST-0098597). 



The Hanle effect is complicated by the fact that three coordinate systems are involved: star, observer, 
and the local magnetic field orientation. A standard approach for solving for the requisite angles is to employ 
spherical trigonometry. More challenging are the interior angles of the spherical triangles in part because 
of quadrant ambiguity and in part because of potential divide-by-zero problems. An equivalent and more 
straight-forward approach is to employ vector relations, which we describe here. 

The givens include the stellar and observer axes and the magnetic field geometry. These immediately 
provide the scattering direction, and vector magnetic field at any point. The stellar axes also allow one to 
specify the direction of the incident radiation throughout the scattering medium. 

For the various spherical triangles of Fig. 3, the "arc relations" are given by dot products, with 



A. DETERMINATION OF ANGULAR QUANTITIES 
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cos^s = f-B, 
cos^i = z-B, 



(A5) 
(A6) 



where z 



Consider the spherical triangle of Fig. 15, drawn from unit vectors d, b, and c. Let the interior angles be 
a, (3, and 7, respectively. Moreover, let the arcs be specified by 6^^, ^ab) and ^bc- Then one can determine 
the cosine and sine of an interior angle by deriving the unit tangent vectors at the vertices of interest. 

For example to determine the angle a, one defines the unit tangent vectors tab and fac at the vertex of 
a in the directions of b and c. Then one has that 
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For the Hanle effect with incident unpolarized radiation, the two interior angles of relevance are ig and 
6. For the first of these, and working through the vector relations, one obtains 
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and for the second angle, one has 
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B. WCFIELDS FOR SLOW ROTATORS 

In Wind Compression theory (Bjorkman & Cassinelli 1993; Ignace et al. 1996), the distortion of the wind 
flow is described kinematically. Ignoring pressure gradient terms, and assuming only radial forces, one can 
derive an expression for streamline flow. The general properties of wind compression have been confirmed 
with hydrodynamic simulations by Owocki, Cranmer, & Blondin (1994), although the effects appear to be 
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inhibited when non-radial accelerations that arise in the line-driving of hot star winds are included (Owocki, 

Cranmer, & Gayley 1996). 

However, retaining the basic elements of Wind Compression theory, Ignace et al. (1998) extended the 
method to allow for "weak" magnetic fields. The key assumptions are that the magnetic field is frozen-in 
and dominated by the hydrodynamic flow. Consequently, the known flow geometry determines the magnetic 
fleld topology. 

Referring to Ignace et al. (1998) for the derivation, the key expression that determines the vector 
magnetic field throughout the flow is 



where 5* is the surface field strength, dji/dno is the "compression factor" that describes how neighboring 
streamlines evolve throughout the wind flow, V is the vector velocity in the co-rotating frame, and is the 
radial velocity. All of the velocity components and the compression factor are knowns. The scalar strength 
of the field is 
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where z;rot is the equatorial rotation speed of the star, and </>' is a "deflection" angle that determines the 
streamline flow and is described in Bjorkman & Cassinell (1993). The zero subscript indicates a value at the 
base of the wind. In the upper hemisphere, streamlines move toward the equator with i? > t^o for r > ii* . 
In the lower hemisphere, streamlines also move toward the equator, now with d < i9q. The field is taken 
to be radial at the wind base (split monopole), but tends to toroidal configuration at large radius, with 
B oc sini?, like the magnetic field that is dragged out in the solar wind. 

The vector orientation of the magnetic field at any point is given by the angles -dB and (fB defined by 
the relations 



Br 



B, 



COSIER = Zi, ■ B = C0S1? sin 



B 



B 



(B3) 



and 



^ By, Br sin -dsmip + B^ cos ip + B-^ cos d sin (p 

Bx, Br sin 1? cos (p — B^ sin ip + B^ cos cos (p ' 

So far, these expressions are valid for all rotation speeds, for which the streamlines are not equator- 
crossing. The assumption of slow rotation with only mild distortions of the wind flow from spherical implies 
the following simplifications (Ignace 1996; Bjorkman & Ignace in prep). 

First, the terminal speed depends only weakly on the initial latitude of a streamline and may thus be 
assumed a constant. 

Second, the streamline parameter becomes (/)' « (^gq(r) sin 6*0, a separable function of radius and initial 
stellar colatitude, with 
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KM = — 



1/7 



5^.(1/7,1-1/7), 



(B5) 



where By is the incomplete Beta function, with Y — 1 ^ vo/vr- 

Third, the radial and latitudinal dependence of the compression factor can be fit with the following 



form, 



1^ ) « (1 + + [(cos^;,)-^ - (1 + ^'.l)-'] sin^ t?. 



where the exponent factor q is 



(B6) 



g(r) =3tan^!)^q(r). 

Consequently, the wind density as a function of radius and latitude is described by 



(B7) 



p{r, jj) = psph 



djjL 



(B8) 



where psph = po •w~^. The upper panel of Fig. 16 shows the variation of the wind density with co-latitude 
at different values of cj)'^^. The lower panel displays the asymptotic equator-to-pole density contrast Peq/Ppoie 
as a function of the ratio Wrot/^'oo- 

Finally, with regard to the magnetic field, it is adequate to assume that i?o ~ which yields for the 
total field strength 



cos^ 1? sin^ (j)' + sin^ i? ( — - ^ 

it* 



1/2 



In the special case that 7=1, the function ^L(r) takes on an especially simple form, with 



(B9) 



(BIO) 



Since 7 = 1 is the wind velocity used throughout our model line profile calculations, Eq. (BIO) used in 
conjunction with the preceding expressions define the flow geometry and the magnetic fleld throughout the 
wind. 
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Fig. 1. — A plot of the Hanle field strength -Bnan for 
resonance scattering lines common to astrophysics. 

Note that a few of these lines arc singlets, but most 
are Li-like doublets. Of the doublets, only the shorter 
wavelength component will show a Hanle effect, since 
the longer wavelength component scatters isotropi- 
cally and produces no line polarization. Evident is 
that shorter wavelength lines, having higher A- values, 
tend to have higher Hanle field values. The shorter 
wavelength lines are also associated with more highly 
ionized atoms that are commonly observed in hot star 
winds. 




% 9 

Fig. 2. — The stellar coordinate system, both Carte- 
sian and spherical. The Z*-axis is taken as the refer- 
ence direction for the field geometry. If the direction 

f corresponds to the illumination ray of a scattering 
point at the dot, then 1} = '&i and (p = ip^. 



Fig. 3. — The scattering geometry centered at the 
scattering point of Fig. 2. The unit vectors r and s 

are for the radial direction and the scattering direc- 
tion. The angles (x, a) are polar spherical angles for 
the observer (along s). The angles 'ds and ^ps spec- 
ify the vector orientation of the local magnetic field, 
B. The angles (6', are polar spherical angles for a 
system defined by this local magnetic field. The an- 
gle is is a Mueller rotation angle that rotates the Q 
and U polarizations from the reference frame of the 
magnetic field to that of the observer's measurement 
axes. 
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Fig. 4. — A sequence of P Cygni profiles produced by 
our code using the escape probability method with 
volume emissivity. The profiles are for optical depth 
scales T; as indicated, with thicker lines producing 
deeper blueshifted absorption and stronger redshifted 
emission. The dip that appears just blueward of line 
center in these profiles is real. In each case the wind 
velocity law is with 7=1, and for an ion species of 
constant ionization fraction. 
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Fig. 5. — Variation of Sq and Su across an expanding 
and optically thin spherical shell that is threaded by 
an axial magnetic field of constant strength. The 
Stokes source functions arc plotted against observed 
Doppler shift, from the front side at {voha/vmax = 
— 1) to the rear side at (wobs/i'max = +!)• The two 
panels at left show the effect of varying the Hanlc 
effect via the parameter cosa2 = 1-0, 0.75, 0.5, 0.25, 
and 0.0 for B/Br = 0.0, 0.9, 1.7, 3.9, and oo. The 
panels at right show the effect of varying the viewing 
inclination (with the solid line for i = 90° and the 
dotted short-dashed line for i = 0°). Note that the 
upper right panel uses a different scale from the other 
panels. 
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Fig. 6. — Results for the Q and U polarized flux pro- 
files from our numerical code for an axial magnetic 
field of constant strength with r; = 1.0 and 7 = 1. 

The sequence of line types arc the same as in Fig. 5, 
but the values used for cosq;2 have been changed 
slightly so that the profiles may more easily be seen. 
The left panels arc for cosz = 0.5 with different a2 
values, and the right panels are for cosq;2 = 0.3, but 
for different values of i. The vertical scale is different 
for the upper right panel. The integrated polarized 
line fiux in both Q and U are generally non-zero. 
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Fig. 7. — Shown as dotted lines are isovelocity zones, 
with the observer to the right. The soUd hnes are 
curves demarcating regions of the wind where the 
Sobolev optical depth in the line is thick (rs > 1 
interior to the curves) and thin {ts < 1 exterior to 
the curves). Three cases are shown, with t; = 0.5, 1, 
and 2. The ti =1 case is shown as a bold line, since it 
is the value adopted for most of our line calculations. 
In this particular case, some portion of every isove- 
locity zone is optically thick in the line. In the single 
scattering approximation, polarized line emission is 
produced only at locations where rs < 1. 
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Fig. 8. — The percent polarization of the Hanlc effect 
for a dipole magnetic field in a spherical stellar wind, 
(a) Upper left: Shown are polarized Stokes Q profiles 
as normalized to the continuum emission outside the 
line frequencies. For a dipole field as seen from the 
side in a line with ti = 1.0, larger Hanle ratios (as la- 
belled) yield higher line polarizations, (b) Lower left: 
As in (a), except now r; is varied as indicated for a 
fixed Hanle ratio. Note the line asymmetry for the 
two lower optical depth cases is a consequence of stel- 
lar occultation. The strongest polarization results for 
fairly optically thin lines, that sample stronger fields 
deeper in the wind, (c) Upper right: Plotted is the 
total polarization P = 



u/rtot in percent. 
The total flux Ftot includes the P Cygni line, which 
is inherently asymmetric in its profile shape. Over- 
all, the effect of viewing inclination is to reduce P for 
perspectives that are increasingly pole-on. (d) Lower 
right: Same line types as in (c) but now for the polar- 
ization position angle Op . Owing to symmetry, U = Q 
at all points in the line for i = 0° and i = 90°. A 
left-right symmetric {/-profile remains for intermedi- 
ate viewing inclinations (which spatially corresponds 
to fore-aft symmetry about the star). 
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Fig. 9. Polarized line profiles for slow magnetic 
rotators at three different values of Wrot/^'oo as indi- 
cated. The magnetic geometries are determined by 
the WCFields model. The wind density and flow ve- 
locity are taken to be spherical, so that the polar- 
ization arises only from the Hanle effect itself. The 
different curves in each panel are for Hanle ratios of 
i3,/i?Han = 0; 3, 10, 30, and 100, with stronger po- 
larizations resulting for larger Hanle ratios. The in- 
clination is 90° and ti = 1. A weak U signal exists, 
but we choose not to show it here. 



WCFields with v„i / v.= 0.08 and -i5= 45" 




Fig. 10. — Using approximations that apply for slow 
stellar rotations, this figure shows the variation of 
the magnetic field components from the WCFields 
model as a function of radius. The individual compo- 
nents are normalized to the local total magnetic field 
strength. These curves are for the case t^rot/i'oo = 
0.08 at a latitude of ?9 = 45° so as to roughly max- 
imize the contribution of the latitudinal component 
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WCFields with WCZ Density 



0.2 


0.2 
^"-0.4 
-0.6 
-0.8 



11 ' ' ' ' 1 ' ' ' ' 1 ' 

: \\ V 30" 


"Tip- 

7 ^ 


r 90' 




l! 


\ - 


: rrr. 
■ 11 \ / \ 


~':\ : 
Ak , - 


: YJ; - 

T , , . . 1 , , , . 1 . 


. , 1 , . , . r 


-1 -0.5 


0.5 1 




Fig. 11. — A series of polarized line profiles for a 
WCFields model with Vrot = 0.08t;oo to highlight the 
influence of optical depth and viewing inclination ef- 
fects. The surface field strength is _B, = lOBnan for 
every profile. The uppermost plot shows the con- 
tinuum normalized Q flux for i = 90° with the line 
optical depth varied as indicated. The four panels 
below show the influence of viewing inclination for 
Ti = 1.0: At left are plots of Fq/F^ and Fu/F^ in 
percent polarization and at right are plots of P (this 
being normalized to the full P Cygni line profile and 
not just the continuum level) and Op. The differ- 
ent curves are for the inclination values shown in the 
panel for Fq. 




Fig. 12.— The Hanle effect for a WCFields model 
with v-cot /voo = 0.08 but now with an asphcrical wind. 
The bold solid line is the line polarization that results 
when there is no magnetic field. The dotted lines 
arc with the Hanle effect (to be compared with the 
middle panel of Fig. 9. The view is edge-on and t; = 1 
for all of the profiles. For B = Q, peak polarization 
occurs near line center and is positive because there 
are relatively more line scatterers in the vicinity of 
the equator, than near the poles. 
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Fig. 13. — Continuum polarizations arising from 
Thomson scattering of starlight in the mildly dis- 
torted WCZ models. The abscissa is v^othao relevant 
for the models presented in this paper. Also shown at 
top is the corresponding asymptotic equator-to-pole 
density contrast. The electron scattering is assumed 
to be optically thin, with envelope optical depths rg 
as indicated. These tq values are for an equivalent 
spherical envelope. The continuum polarizations are 
for an edge-on view, and will scale as sin^ i for other 
viewing inclinations. The filled dots indicate the cor- 
responding models used in the line calculations for 
the WCFields magnetic geometry assuming a spher- 
ical wind density. 



WCFields with Initial Magnetic Monopole 
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Fig. 14. — Polarized line profiles in Q and U for a 
WCFields model, with Vrot/vao = 0.08 and a spheri- 
cal density, but now with an initial magnetic field at 

the wind base that is a magnetic monopole. Com- 
pared to the basal split monopole configuration used 
in Fig. 9, the Q polarization has dropped, and a sig- 
nificant and antisymmetric U profile is now evident. 
Consequently, the U profile, when the polarization 
measurement axis is aligned with the source symme- 
try axis, is sensitive to the symmetry of the circum- 
stellar magnetic geometry. 
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Fig. 15. — Specification of spherical triangle unit vec- 
tors, spherical polar angles, and interior angles used 
in relating the evaluation of the angles by vector prod- 
uct relations. 
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Fig. 16. — The wind geometry of WCZ models using 
our approximations for slowly rotating stars. The 
upper panel shows the variation of the normalized 
wind density with co-latitude at different values of 
(/>gq. Below is the asymptotic wind density contrast 
between the equator and the pole as a function of 
the stellar rotation normalized to the wind terminal 
speed. 



